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REVIEW  Ur  INVEST  i  GAT  I O::  PROGRESS 


A .  Ini t ial-Va la e  Problems  in  Stratified  Si lear  Flow 

The  question  of  the  stability  of  a  stratified  shear  flow  is  addressed 
through  the  investigation  of  the  initial-value  problem  defined  for  a  two-layer 
fluid  of  infinite  extent  with  uniform  velocity  and  density  in  each  layer. 
Solutions  of  this  problem  have  application  to  the  generations  of  unstable 
motions  in  the  atmosphere. 

The  eigenvalues  of  this  problem  can  be  expressed  as  the  solutions  oi  a 
quadratic  polynomial  in  the  frequency  with  coefficients  which  are  functions  of 
the  wavenumbers,  and  relative  densities  and  velocities.  The  value  of  the 
discriminant  of  this  equation  determines  if  the  system  is  stable  or  unstable. 
Figure  1  is  the  locus  of  the  neutral  values  of  the  discriminant  as  a  function 
of  velocity  and  wavenumber.  The  graph  indicates  that  there  is  a  threshold 
velocity  below  which  the  system  is  stable.  Figure  2  is  an  example  of  the  values 
of  the  imaginary  part  of  the  complex,  frequency  for  waves  traveling  at  angles 
of  obliquity  of  0°,  30°  and  60°  to  the  mean  velocity. 

Solutions  to  the  initial-value  problem  driven  by  a  Gaussian  pulse  showing 
the  initial  distortion  and  tiie  asymptotic  development  are  shown  in  Figures  3 
and  4,  respectively. 

The  initial  distortion  is  evaluated  using  a  power  series  expansion  and 
approximation  of  the  complex  dispersion  relation  as  a  polynomial  in  the  hori¬ 
zontal  wavenumbers.  The  expansion  for  non-dimensional  Lime  parameter  0,  0.03 
and  0.75  are  shown  in  Figure  3. 

The  two-dimensional  asymptotic  expansion  using  saddle-point  techniques 
for  this  non-conservative  system  is  evaluated  at  values  of  the  time  parameter 
of  20  and  60  and  is  shown  in  Figure  4.  It  is  apparent  that  the  system  tends 
to  become  two-dimensional  (long  standing  wave  patterns)  in  the  large. 


A  complete  report  el  this  work  is  contained  in  the  Ph.D.  dissertation 


of  J.E.  Bradt  with  a  version  for  publication  to  be  submitted. 

A  second  piece  of  research  as  been  completed  (submitted  for  publication) 
and  concerns  the  question  of  other  singularities  in  the  eigenvalue  spectrum. 

A  recent  paper  by  Chimonas  (G.  Chimonas,  1970,  Journal  of  Fluid  .Mechanics ,  90. . 
1-69)  drew  attention  to  the  fact  that  essential  singularities  are  present  and 
thereby  result  in  (a)  algebraic  growth  in  time  for  the  perturbations  and  (b) 
are  not  form  preserving.  Unfortunately,  Chimonas  did  a  purely  inviscid  calcu¬ 
lation  and  our  r.n.nlyri-  a'.icva  ‘-u-t  over,  a  small  bit  cf  visceaii,  (large  Reynold 
number)  negates  this  result.  Nevertheless,  the  general  Laplace  transform 
technique  does  reveal  that  considerable  care  must  be  given  to  perturbation 
problems  in  stratified  shear  flow. 

The  general  study  of  the  problem  of  shear  flow-internal  wave  interaction 
continues  as  well.  A  modelled  situation  is  proving  quite  productive,  partic¬ 
ularly  when  the  major  goaL  is  a  representative  of  a  spectrum  of  the  fluctuating 
energy.  Specifically,  tne  work  includes  (a)  viscous  effects,  (b)  a  constant 
mean  shear,  (c)  a  constant  Brunt-Vasailia  frequency  (exponential  density 
stratification),  but  omits  any  effects  due  to  boundaries.  This  last  point 
seems  reasonable  enough  in  view  of  the  fact  that  we  are  concerned  with  motions 
in  the  interior  of  the  atmosphere.  .'O--  ver,  the  technique  used  for  solution 
is  novel  and  is  capable  of  producing  results  for  arbitrary  three-dimensional 
disturbances  with  the  aid  of  an  analogue  computer.  Both  a  Master's  thesis  and 
a  publishable  version  of  this  work  will  soon  be  available. 


B.  Fluctuations  i_n  Geophysical  Boundary  Layers 

Concentration  on  the  analysis  of  perturbations  in  a  turbulent  Ekman  layer 
has  been  the  center  of  this  research  since  the  last  report.  Besides  desiring 


the  characteristics  of  the  structure  in  the  physical  situation,  an  answer  to 


the  important  question  of  body  force  effects  on  turbulence  was  sought.  In 

order  to  reduce  the  system  to  a  workable  (and  understandable)  means  it  was 

decided  to  allow  the  important  parameter  to  be  the  relative  scale  depth  of  the 

logarithmic  portion  of  the  boundary  layer  (near  wall)  to  that  of  the  total 

shear  zone  (including  spiral).  When  this  ratio  is  too  small,  the  layer  is 

unstable,  a  result  quite  different  from  conventional  flat  plate  turbulent 
i 

boundary  layers  where  it  is  v;ell  known  that  a  fully  turbulent  layer  is  stable. 

The  ramifications  of  this  are  being  explored  further. 

A  report  on  both  the  laminar  (wave  packets)  and  turbulent  Ekman  layers 
will  be  available  with  the  Ph.D  dissertation  of  r..F.  Spooner. 

C .  Internal  lv_  or  Ext  email  y  _Pe  rturbed  Boundary  Layers  over  a^Flat  Flare 

We  consider  steady  viscous  incompressible  flow  over  a  semi- inf inite  flat 
plate  under  the  influence  of  either  internal  or  external  time-dependent  per¬ 
turbations.  Experimentally,  an  internal  perturbation  may  be  realized  by  a 
vibrating  ribbon  and  an  externa)  perturbation  corresponds  to  small  fluctuations 
superimposed  on  a  uniform  frcestream. 

We  show  that  in  either  case,  when  the  amplitude  of  perturbations  o  is  less 
than  the  boundary  laver  thickness  z  but  larger  than  e“  the  unsteady  correction 
terms  to  the  conventional  boundary  layer  equations  are  governed  by  certain 
linear  variational  equations  with  variable  coefficients;  these  coefficients 
involve  the  Blasius  velocity  functions.  For  the  case  of  internal  perturbations 
we  obtain  certain  unsteady  boundary  conditions  to  be  imposed  on  the  flat  plate 
while  for  external  perturbations  the  unsteady  boundary  conditions  are  specified 
at  infinity  in  the  vertical  direction. 

The  derivation  of  these  results  proceeds  systematically  from  the  N’avier-Stokes 
equations  using  conventional  asymptotic  arguments  and  does  not  involve  the  ad  hoc 


assumption  of  parallel  flow  in  the  boundary  layer.  In  fact;  both  the  horizontal 
and  vertical  velocity  components  of  the  Blasius  solution  occur  in  our  variational 
equa  tions . 

A  preliminary  program  for  numerically  solving  the  variational  equations 
for  the  case  of  a  sinusoidal  internal  perturbation  is  complete,  the  results 
are  qualitatively  in  accordance  with  experimental  observations  in  that  the  dis¬ 
turbance  region  diffuses  in  the  vertical  direction  and  is  convected  downstream. 
This  numerical  procedure  is  only  valid  for  moderate  frequencies  and  does  nor 
apply  for  the  more  interesting  case  of  high  frequency  disturbances.  To  analyze 
this  limit  we  have  rescaled  the  variables  in  our  variational  equations  and  we 
are  currently  working  on  a  second  numerical  code. 

Paralleling  these  numerical  studies  we  have  derived  an  asymptotic  approxima¬ 
tion  of  the  perturbation  equations  in  the  limit  of  large  excitation  frequency. 

In  tills  limit  the  basic  problem  reduces  to  a  weakly  nonlinear  diffusion  equation 
with  an  oscillating  boundary  condition.  An  approximation  valid  for  short  times 
lias  been  derived  anu  indicates  that  the  disturbance  in  the  boundary  layer  grows 
with  time.  In  order  to  derive  a  uniformly  valid  result  for  long  times  we  plan 
to  develop  a  multiple  variable  solution.  Multiple  variable  expansions  have  been 
used  with  much  success  in  wave  propagation  problems  as  well  as  in  some  examples 
of  weakly  nonlinear  elliptic  equations.  However  there  are  no  significant  results 
available  for  weakly  nonlinear  parabolic  problems  with  time  dependent  boundary 
conditions.  We  propose  to  concentrate  on  this  point  during  the  next  grant 
period . 

An  initial  report  on  this  work  will  be  given  at  the  forthcoming  meeting  of 
the  Division  of  Fluid  Dynamics,  American  Physical  Society,  to  be  held  at  the 
University  of  Notre  Dame  in  November  1979. 

The  alternative  approach  to  analysis  of  the  continuous  spectrum-earlv  time 


problem  for  the  boundary  layer  is  stiLl  is  progress. 


Basic  equations  are 


available  with  a  means  for  solution  still  not  complete.  This  work  is  of  high 
priority  . 

D.  Report  on  a  Turbulence  Initiation  (.aicula t iona 1  Procedure 

Analytical  and  experimental  results  obtained  by  a  number  of  investigators 

(see,  for  example,  Landahl,  M.T.,  J FM  j>6,  1972,  p.  775;  Landahl,  M.T.  and 
i 

Criminale,  W.O.,  J FH  79  ,  1977  ;  Kvgnanski,  I.,  Haritonidis,  J.H.,  and  Kaplan,  R.E., 
JFM  92,  1979,  p.  503;  Gaster,  M.  and  Grant,  I.,  Proc.  Roy.  Sue,  bond .  3^*7 ,  1975, 
p.  253;  Kim,  H.T.,  Kline,  S.J.,  and  Reynolds,  l.’.C.,  J F'?  50 ,  1971,  p.  133,  and 
the  reterences  contained  therein)  help  focus'  attention  on  both  the  importance 
ana  the  complexity  of  the  nonlinear  processes  which  take  fhe  place  after  the 
first  appearance  of  a  turbulent  spot  in  a  boundary  layer.  It  appears  useful  to 
attempt  to  .supplement  those  investigations  with  a  reasonably  accurate  computa¬ 
tional  approach  to  the  problem,  and  tor  this  purpose  we  have  considered  several 
possible  techniques,  and  have  experimented  briefly  with  a  program  based  on  one 
of  them.  The  program  appears  promising,  but  it  is  clear  that  much  remains  to 
be  done. 

Using  Cartesian  index  notation  (and  the  usual  summation  convention),  the 
Navier-Stok.es  equations  for  incompressible  flow  may  be  written 
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wiiere  iu  is  the  k1  tli  component  of  velocitv,  p  is  the  density,  ",  the  component 
K  K. 

of  body  force,  and  v  the  kinematic  viscosity.  Taking  the  curl  of  Eq .  (1)  we 
obtain  after  simple  manipulation  an  equation  for  the  vorticitv  vector  : 
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and  we  note  also  that 
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The  advantage  of  using  a  vorticity  equation  is  that  the  pressure  derivative 
is  eliminated.  Our  reason  for  including  a  body  force  term  is  twofold.  First, 
its  use  provides  a  convenient  and  flexible  triggering  mechanism,  which  avoids 

I 

the  necessity  of  inserting  a  perturbation  in  either  the  entering  or  existing 
laminar  flow.  Second,  there  is  experimental  evidence  to  the  effect  that  a 
longitudinal  flow  vortex  can  play  an  important  role  in  the  growth  of  turbulence 
near  an  initial  turbulent  spot:  one  way  in  which  such  a  vortex  can  be  generated 
is  if  there  is  an  unequal  division  of  flow  around  a  turbulent  spot  (the  instab¬ 
ility  of  parallel  channel  flows  is  well  known),  and  this  kind  of  inequality 
can  be  easily  triggered  by  the  use  of  an  appropriate  small  body  force  distribu¬ 
tion. 

Our  approach  is  to  take  time  steps  in  vorticity,  using  Eq .  (2),  and  to  then 
use  Eq .  (3)  to  determine  the  new  velocity  field.  The  region  of  interest  is  any 
simple  geometrical  region  containing  the  turbulent  spot,  with  the  boundaries 
sufficiently  far  removed  that  in  the  initial  portion  of  the  turbulence  build-up 
there  will  only  be  a  negligible  perturbation  felt  at  these  boundaries.  This 
will  of  course  not  be  precisely  true,  and  in  fact  it  will  not  be  at  all  correct 
on,  say,  a  flat  plate  boundary  near  which  the  turbulence  is  being  generated. 
Consequently,  a  correction  procedure,  near  anv  such  boundary,  will  be  necessary. 

To  take  a  time  step  in  Eq .  (2),  we  have  so  far  used  a  simple  explicit 
procedure.  Although  this  can  be  chosen  to  be  both  stable  and  accurate,  and  has 
in  fact  been  quite  satisfactory  for  the  initial  ’ nves t i ga t ion ,  it  is  recognized 
that  larger  time  steps  are  feasible  for  implicit  procedures  of  the  kinds 


developed  for  diffusion  type  equations.  Possible  improvement  in  this  particular 


procedure  will  be  examined.  Using  the  new  vorticity  values,  Eq  .  (3)  has  the 
form  of  a  Poisson  equation  for  each  component  of  velocity.  Assuming  first 
that  the  boundary  values  of  velocity  are  those  of  the  laminar  flow,  new  interior 
values  of  velocity  corresponding  to  the  new  vjrticities  nav  be  obtained  by  •- 
successive  optimal  over-relaxation  procedure  which  (for  a  prismatical  region, 
say)  is  very  efficient.  The  number  of  iteration  sweeps  required  to  reduce  the 
error  by  a  factor  of  100  is  of  the  order  of  the  number  of  mesh  points  alone  a 
side;  the  initial  error  is  kept  small  by  using,  for  a  first  guess,  either  the 
velocity  values  appropriate  to  the  previous  time  step  or  extrapolations  from 
those  values . 

The  coupling  between  Eqs.  (2)  and  (3)  required  for  boundary  corrections 
arises  at  the  next  step  in  the  program.  In  order  to  compute  interior  vorticity 
values  for  the  next  time  step,  we  need  new  vorticity  values  on  the  boundarv. 
However,  we  do  not  have  a  priori  values  for  velocity  derivatives  (.and  hence 
vorticities)  on  the  boundary;  thus  we  must  extrapolate  to  obtain  these  values 
from  interior  values.  This  kind  of  calculation  is  known  to  lead  to  possible 
inaccuracies  and  even  instabilities,  and  it  may  be  necessary  to  devise  safe¬ 
guards  against  such  inaccuracies.  There  are  some  situations  in  which  analvtical 
inturmatiuii  n.vuiviin;  ucu.ativc  •<#«*.'  are  availab2°  from  the  above  equations, 
and  if  so  the  numerical  calculations  can  be  simplified. 

Another  problem  which  must  ultimately  be  dealt  with  is  that  of  mesh  refine¬ 
ment  near  critical  regions,  for  the  purpeae  cf  '’ah  in''®'1  '.r—rir-,  Our  experience 
with  other  problems  of  diffusion  or  potential  flow  character  leads  us  to  expect 
this  possibility  to  be  useful  (in  that  it  reduces  the  amount  of  computer  time 
required  for  a  desired  accuracy  of  result),  we  have  not  experimented  with  this 
feature  in  the  present  problem. 

A  preliminary  program  has  been  written,  and  applied  to  the  initial  phases 
of  turbulence  growth  near  a  flat  plate  boundary.  The  mean  stream  flow  here  was 


Poiseuille  flow  character,  with  an  initially  parabolic  veiocitv  distribution. 
Ncn-diraenpi onaiizati on  was  used,  and,  in  terms  of  the  distance  to  the  inflection 
point  in  the  parabolic  profile,  the  Reynolds  numbers  used  were  primarily  in  the 
range  of  a  few  hundreds  to  a  few  thousands.  Disturbances  were  generated  by 
impulsive  localized  body  force  contributions,  and  with  about  L50u  mesh  points, 
a  few  tens  of  time  steps  were  taken. 

With  the  limited  time  and  computational  facilities  available,  onlv  verv 
i 

preliminary  resuLts  could  be  obtained.  However,  the  program  did  appear  to  be 
both  accurate  and  stable,  at  least  for  the  initial  part  of  the  motion.  To  date, 
both  the  space  and  time  meshes  are  crude,  the  control  of  perturbed  boundary 
values  is  primitive,  and  —  as  previously  mentioned  —  refinements  aimed  at 
increasing  accuracy  or  stability  have  not  been  made.  Our  general  feeling  is 
that  the  method  appears  to  be  well  worth  further  investigation  as  a  tool  in 
better  understanding  the  growth  of  turbulent  regions,  and  we  propose  to  under¬ 
take  such  an  investigation. 
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Initial  distortion  of  Caussian  pulse  of  standard  deviation 


